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' Abstract 

o 

This paper shows that in intensity correlation measurements there will be clear and un- 
ambiguous signals that new-physics particles are, or aren't, parabosons. For a parabosonic 
field in a dominant single-mode, there is a diagonal P p -representation in the \a evenoc id > 
coherent state basis. It is used to analyze zero-time-interval intensity correlations of para- 
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bosons in a maximum-entropic state. As the mean number of parabosons decreases, there 
is a monotonic reduction to - of the constant bosonic "factor of two" proportionality of the 
second-order versus the squared first-order intensity correlation function. 
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Statistics based correlations, in particular second and higher-order intensity correlations [1,2], 
have proven to be very useful in interferometry measurements and in dynamical investigations of 
new phenomena in areas ranging from astrophysics to nuclear and high energy physics [3]. It has 
been shown that relativistic, local quantum field theory allows particles which obey parastatistics, 
an elegant generalization of bose and fermi statistics based on the permutation group [4] . While 
a use for parastatistics in the description of natural phenomena has yet to be found, a wealth of 
new but poorly studied discoveries have been recently made in astrophysics/cosmology including 
an inflation era, dark matter, and dark energy. New particles are expected to be produced in 
experiments at the Large Hadron Collider such as Higgs particles associated with electroweak 
symmetry breaking (EWSB) and supersymmetry particles associated with the energy-scale in- 
equality A EWS b << ^Planck- Therefore, for incisive analyses, it is both timely and important to 
know both the selection rules and the novel statistical effects/correlations [4,5] for the purpose 
of detecting the presence of effective or underlying fundamental quanta obeying parastatistics. 
This paper shows that by intensity correlation measurements there will be definitive signals that 
new-physics particles are, or aren't, parabosons. 

For single- mode parabosons (pB's), we introduce an index i = even, odd . The single- mode pB 
commutation relations [4,6] are [a, {a*, a}] = 2a and its adjoint, with N B = H B — | = \{a) , a} — |. 
In the occupation number basis, the states are 
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with n e = 2N, n Q = 2N + l;iV = 0,l,...; and where the p-factorials are defined 



(n e ) p ! = (2N) P \ = 2N(2N -2+p){2N -2)---4(2+p)2p; (0) p ! = 1 
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Note (0) p ! = 1, (l) p ! = p, and (n ) p \ = (2N + l) p \ = (2N + p)(2N) p \. Mathematically, in counting 
the n p integers, to reach the next odd (even) n p integer, one adds p(l) to the usual integer. For 
fixed-order p = 1,2,... parabosons, this {n) p \ product of n factors is analogous to the Gibbs- 
factorial n! for ordinary bosons (p = 1). For creation of one more pB, the next to vacuum position 
has p(l) openings when n even (n odd ) parabosons are already present; the other positions have 1 
opening. 2 To make an odd number of insertions depends on p but to make an even number 
does not. Because , a are raising and lowering operators, all pB single- mode a\ a commutation 
relations in this \n efi > basis follow from (1). 

For single-mode parafermions (pF's), in the occupation number basis, there is a band of 
states. The pF commutation relations are [c, [c^c]] = 2c and its adjoint, with N F = H F + | = 
\[c\ c] + |. The pF number states are equivalently 



2 This characterization follows from [a, (at) 2 ™] = 2n( a t) 2 ™- 1 and [a, ( a t) 2 ™ +1 ] = (a^) 2n (2n + [a, at]) since 
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with c|0 >f= c^\p >f= 0. Analogous to the (n) p \ for pB's, for pF's 

{n} p \ = {n(p -n+ l)}{(ra - l)(p - n + 2)} • • • {2(p - !)}{lp} 



(5) 



' (p-n)! 

which is a product of n bi-factors. Note {0} p ! = 1 and {1} P ! = p. In constructing a specific \n >f 
state from the adjacent \n — 1 >p state, the additional {n{p — n + 1)} bi-factor is the product of 
the number of ways to insert the extra " " in the first expression in (4), times the number of 
ways to remove a " c " in the second expression. The Hamiltonian defines the direction |0 >^ to 
\p > F . 

There are important consequences of (n) p \ and {n} p \: For a number "n/' of parabosons initially 
in the mode, we let Prob B (^ ~* ^) be the probability that in some fixed time interval a system 
is found to make a transition 3 from a state "a" to a state "6" by the emission or absorption of a 
single pB 7 . For initially n even = 2N pB's, the decay probability is 

whereas for excitation from the state "6" to "a" , in the same time interval, the absorption proba- 
bility is 

Prob B (— -> — - — ) = (2N) A (7) 
K 2N 2N-1 J K J K J 

By time-reversal invariance, the constant A is the same for emission and absorption. Hence, 
versus ordinary bosons, for the system initially in an n e t, en -mode, there is a p-dependent enhanced 
7 stimulated-emission. Likewise, for initially n dd = 2iV + 1 pB's, for emission/absorption 

Prob B ( „^ - -> — h -— ) = (2N + 2) A; Prob B ( n J : -> 7^7) = (27V + p) A (8) 



3 A «. 



A "soft accent" denotes a single para-particle. 



Hence, for ordinary bosons, p — 1, there is stimulated emission. However, for the system initially 
in an n odd -mode, for p = 2 the probabilities are equal; but for p > 2 pB's there is a p-dependent 
enhanced 7 absorption versus emission. In summary, from consideration of single-mode pB 
statistics versus bose statistics, for a system in an initial pB n even - mode ( n^-mode ), there 
is a p-dependent enhanced 7 stimulated-emission by the system (stimulated-absorption by the 
system); otherwise there is no p-dependence in these transition probabilities. This is different 
from the generic bunching signatures of ordinary bosons! 

Also, in contrast, for parafermions the statistics-favored transition probabilities are towards 
the "half-full/half-empty" \n mi d > parafermion-state(s) at the middle of the pF band of "p + 1" 
states, due to the bi-factors and time- reversal invariance. From {n} p \, for p even there is a single 
\n m id >= > state but for p odd there is a pair of mid-band states \n mi ^hi g h/iow = Podd±1 >. 

For ordinary fermions, there are only these two mid-band states: 

For a number "raj" of pF's initially in the mode, we let Probp(^j — > ^) be the probability 
that in some fixed time interval a system is found to make a transition from a state "a" to a state 
"6" by the emission or absorption of a single pF v . If initially rij > is in the less-than-half-filled 
part of the pF band with raj < ( £s ^ 2 -, Pod ^~ 1 ), then the ratio of probabilities for emission of a single 
pF v by the system versus v absorption is greater than one, 

P^M^^^i) nfa-m + l) l 2' 2 ; ^ 

On the other hand, if initially rij < p is in the more-than-half-filled part of the pF band with 
nj > [Sssssl^ Podd+i ^ then the ratio of probabilities for emission of a single pF v by the system 
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versus v absorption is less than one, 



, / Peven Podd 1 \ 

< 1; p > ri! > (— ^- , ) 



(10) 



As the pF-mode approaches the the middle of the pF band, this is monotonically a smaller 
statistics effect. Simultaneously, if one ignores the differences in the A factors, the v transition 
rates themselves increase as the pF-mode approaches mid-band since the mid-band versus end-of- 
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= \{Podd + 2 + -M > 1; p odd > 1 
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The pB coherent states [6] can be written in terms of a p-exponential function 

oo 

e p (x) = E ?fy 

n=0 v n 

= e e (x) + e a (x); e £j0 (-x) = ±e e , (x). 
In terms of the modified Bessel function I v (x), the even and odd terms are 

e e (x) = (f)^r(|)7^(x);e (x) = (f)^r(^)7 f (x) 

For p—1, e e ,o{ x ) coshx, sinhx. 

Thereby, with a complex- valued, a\a >= a\a > for 



with 



= /P e (|a| 2 ) \a e > + V /P c (|a| 2 ) \a Q > 



lae> ^W)^oW !l 



\a > 



00 |„|2JV+1 , 

V — 27V 4- 1 > 



(12) 



(13) 



(14) 



(15) 



< a e \a e >=< a \a >= 1; < a e |o; >= 



In (14), the important n e>0 coherent-state- mode probabilities are 

These range monotonically in \a\ 2 from P efi — > 1 — ^j-, ^- respectively as |a| 2 — >■ 0, to P efi — > | 



for |a| 2 >> 1, 2 



2 

Although a 2 \a efi >= a 2 \a efi >, it is with a different normalization 



\a + >= v/2P e (|a| 2 ) |a c >; |a_ >= ^(M 2 ) |a > (17) 
that a|a ± >= a|a T >. Note that |a± >= ^[|ct > ± | — ct >], 

| - a >= JP e (\a\ 2 ) \a e > -JP (\a\ 2 ) \a a >, and \a e>0 >= 1 [\a > ± | - a >]. 
There is an associated p-Poisson distribution function, x = \a\ 2 , 

x ni 1 

V„(ni, x) = - — — - — — = I < a > I 2 ; i = even, odd (18) 
(7ii) p ! e p (a:) 

for the probability to be in the nth number-state if the system is in the coherent state \a >. 

V p (n i ,x)^^{l-\f + 0(\a\*)}, 

- 2j W ^ 2( Zl? ]) *m-w\ 2 w + oc^)}, h 2 » i,i 

(19) 

The corresponding p-Gaussian ( p-normal ) distribution approximation to (18) is given in the 
appendix. The p-dependence in the p-Gaussian distribution for the pB coherent state arises only 
through the mean \i and the standard deviation a. This is unlike the explicit p-dependence, and 
the explicit n even versus n dd differences, in the above p-Poisson distribution. From (18), V p 2 ^^ = 
2N +p an d V p 2 (2N)^ = This means f° r the distribution of \n >'s in the coherent state \a > that 
versus the probability of an arbitrary \n even > number-state, the next "odd" number-state is less 
probable for p > 1 than for p = 1, but the probability for the previous "odd" is not p-dependent. 



For two arbitrary coherent states, 



< a\(3 >= ^/P e (\a\*)P e (\m < a e \(3 e > + ^fp^^WW) < » \Po > (20) 



with 



< a e \(5 e >= 6e{a * P) ■ < a \/3 >= 6o{a * P) ■ < a e \(3 >= (21) 
So for large arguments: >> 1, 

| < a|/J > H exp(-i|a - + 0( ^, ^L, ^)} (22) 

and for small arguments: |a|, « 1, 

| < a|)9 > | - {1 - + 0(|a| 4 , |/3| 4 , |a| 2 |/5| 2 )} (23) 

As a consequence of the completeness relation 

I = — I d 2 a fi e (\a\ 2 ) \a e >< a e \ H — / <i 2 a yU (|a| 2 ) \a >< a Q \ , (24) 

TV J IT J 

where in terms both types of modified Bessel functions 

/i e (|a| 2 ) = \a\ 2 K v _^(\a\ 2 )h_^(\a\ 2 )] /i (M 2 ) = |a| 2 ir f (|a| 2 )/ f (|a| 2 ), (25) 

and of the \a > and \f3 > non-orthogonality, the pB coherent states are linearly dependent. They 
are overcomplete. This relation 4 (24) follows by using 



r e (n e + 1) = (2N) P \ = |!|/-rf(|a|2)ir £ _2(|a| 2 )|a| 2 (i + -) 
F (n + 1) = (27V + l) p ! = / °°d(|a| 2 )^ $ (|a| 2 )|a| 2 <«+»-> 



(26) 



4 By substitution of (17) into (24), eq(2.63) of 3rd paper in [6] is obtained. 



This analytic integral representation for the two p-factorials generalizes Euler's formula for T(x). 

oo 

Uses for generalizations of other functions of integers, e.g. ( p (k) = J2 ( n p)~ fc 5 k > 1 , remain to 

n p =p 

found. 

In this \a ei0 > coherent state basis, there is a diagonal P p - representation 5 for the density 
operator p describing the state of the system 

1 I' 1 I' 

p = - d 2 a /i e (|a| 2 ) \a e >< a e \ $ (e) (a) + - / d 2 a /i G (|a| 2 ) \a a >< a Q \ $ (o) (a) (27) 

For p = p , & e, °\a) are real. The normalization condition from p = p e + p Q is 

Trp efi = J e>0 = - fd 2 a /i e , (M 2 ) $ (e ' o) (a) (28) 

7T J 

For instance, if the system is in the coherent state \f3 e > , then pp e = \f3 e >< f3 e \ for & e \a) = 
j^r^ia - (3), $(°)(a) = where 5 2 (a - (3) = 5(Re[a - (3})5{Im[a - (3\). 

To describe a field theoretic system in a maximum-entropic state, we proceed as in the treat- 
ment of ordinary bosons [2] for which such a field-state is often called a chaotic or thermal state: 
In terms of the paraboson number operator N = N B , we maximize the entropy S = — kTr(phip) 
to obtain 

1 ( <N> \ N 

PmaxS = ~~ ~ ~~ 29 

1+ < N > \1+ <N >J 
where < iV > is the mean number of parabosons in the maximum-entropic state. 

Defining r = 1 + I , from the above P p -representation, the corresponding non- negative 

maximum-entropic & e '°\a) functions are 

(e) ri K P _2{r\a\ 2 ) r f Kv_(r\a\ 2 ) 

< N > Kp=i\\oi\ ) <N>J<E{\a\ 2 ) 



5 This is the direct generalization of the bosonic P-representation, sec [7]. We use $( e '°)(a) to denote the weight 
functions, so as to avoid confusion with other functions. 



For ordinary bosons $ [ Z xS (a) = $^ x5 (a) = exp(-|a| 2 / < N >). 



By using these results, in analogy with the scalar-field treatment for ordinary bosons, we 
can analyze Hanbury-Brown and Twiss intensity correlations [1,2] for a paraboson field which is 
in a single-mode maximum-entropic state: We consider zero-time-interval correlations. Using 
PmaxS = PmLs + PmL s > the first-order coherence function is 



G^(0) = Tr[p maxS E(-\ Xl )EM( Xl )} 
= G«(0) + G< 1 >(0) 
= c<N> 



(31) 



p+2<jV> 
l+2<iV> 

where c is a constant factor. The second-order, or intensity, correlation function is 

G( 2 )(0) =Tr[p max5 E(-\x 1 )E^(x 2 )EM(x 2 )EM(x 1 )] 
= 2(c<N>) 2 



p+2<N> 
1+2<N> 



(32) 



If we write a proportionality G (2) (0) = X p [G (1) (0)] 2 , then X p = 2 



!+2<jV> 
p+2<N> 



which shows that as 



the mean number of parabosons decreases, there is a monotonic reduction to - versus the usual 
bosonic "constant factor of two" intensity correlation effect. 



Similarly, we obtain for the next higher-order correlation functions 



G (3) (0) =3!(c<iV>) 3 
G {4) (0) =4!(c< iV>) 4 



ap 2 + bp + c 
.3(1 + 2 < N >) 3 
ap 2 + bp + c 



3(1 + 2 < N >) 3 _ 



(33) 
(34) 



where 



a = l + 2 < N > 

b = 2(l + A<N>+Q<N> 2 ) 

c = 8<N>(l + 3<N>+3<N> 2 ) 



(35) 
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As in (31,32), the last factor in (33,34) depends on < iV >; it approaches "1" for < N > large, 
and a p-dependent value " " as < N 0. 

Formulae for arbitrary n order G e ^ n \0) are listed in the appendix. There is no p-dependence 
in zero-time-interval, intensity correlation functions in the large < N > limit. However, there is 
significant p-dependence as the mean number of parabosons decreases. 

In summary, this paper shows that in intensity correlation measurements there will be clear 
and unambiguous signals that new quanta are, or aren't, parabosons. It will be a complete 
measurement, because the order p of their parastatistics will be determined! 

One of us (CAN) thanks colleagues at Binghamton University and elsewhere, in physics and 
in mathematics, for discussions. This work was partially supported by U.S. Dept. of Energy 
Contract No. DE-FG 02-86ER40291. 

Note added in proof: Recursion relation for intensity correlation functions: 

There is a simple recursion relation between the n even -oidei and the lower- adjacent n 0( M-order 
intensity correlation function 

G (ne) (0) = n e (c<N>) G (ne_1) (0). 

This relation is independent of p and so it might be empirically very useful, for instance for 
measurement of the product of the constant factor and the mean number of parabosons, (c < N >). 
This relation follows from the hypergeometric function expressions in the appendix. 
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Appendix: 



G( n )(0) Formulae for Maximum-Entropic State: 

For G (n) (0) = Gi n) (0) + G^(0) with r = 1 + ^jU, there are integral representations 



oo 



G< n °)(0) = ^fdx x n ° +l KzJrx) Ir(x) 

y ' <N> q 2 ±V/ 2 V/ 

P OO 

G("»)(0) = J dx x n ° +1 Kv_(rx) It^x) 

< N > o 22 (36) 

P oo 

G< ne >(0) = ^Jdx x n * +1 KE-^rx) Iz-!(x) 

V oo 

GK)(0) = J£^jdx x nE+1 Kp(rx) Ie(x) 

<N> g 2 V y 2 V y 

In this appendix, we suppress the overall (c)™ 4 factors. The G^(0) can be written in terms of 
the hypergeometric function or as an infinite series 

r*{n n )(n\ _ . 2 "° 1 ( -2 + ~+2- >1 l -~ + 2^ p _i_ n D _|_ 1 n D _|_ 3. p _|_ -i -2\ 
G e H U J - < ^ >r „ D +2 f(f+l) 2^ll2 + T + 2' T + 2' 2 + X ' r > 

2™° ^ r(|+Af+l+m)r(Af+2+m) _ 2m 



<A r >r n +2 ^ r(|+l+m)r(m+l) 

r'tO/TTl - .. 2 "° r (f + lf + §) r (ir + 5) P i no i 1 n D i l.p -2\ 
G ° V U J - <^ >r n +l r(f) 2^ll2 + — + 2'T + 2' 2' r J 

2"° ^ r(|+AT+l+m)r(jV+l+m) __ 2m 

~~ <Af>r«o+i r(§+m)r(m+l) r 

l U J <7 v >r n e+ i r(f) 2^H 2 ^ 2' 2 ^ 1 >2>' J 

2"° v r(£+jV+m)r(jV+l+m) 2m 

<iv>r« e +i „^ r(§+m)r(m+i) 
G ° W - <7 v >r n e+2 rxf+i) 2^H 2 + — + !, — + !, 2 + 1, r J 

00 r(|+Af+l+m)r(Af+l+m) „ 2m 



(37) 



<Ar>r n e+ 2 ^ Q r(§+l+m)r(m+l) 



For < iV >^ oo, we assume the large x = \a\ 2 region of the integrand dominates, so 
G^}(0) — > |n! < N > n because of the absence of //-dependence in I u (x) and K v (x) for large x. 
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Expansions in p follow from the infinite series expressions, 



^° } (0) = gg*gg£ fr" +1 [l + 2 < N >) 



-p N (N + 1)[N + 2N < TV > +2(/V + 1) < N > 2 ] + • • •} 



G^KO) = 2 -^§^ {p N (N + 1)[2 < N > +2 < N > 2 } + ■ ■ ■} 



+p JV (iV + 1)[7V + 2(iV + 1) < N > +2(N + 2) < N > 2 } + ■ ■ ■} 

(38) 



and 



G (n e )( ) = 2^V!<iV>^+i r ATM + 2 < TV >1 



+p iV - 1 (jV)[7V + 1 + 2(/V + 1) < N > +2(N + 1) < N >< 



Gt^) = ^gg^Kll + 2 < N >] 



+p JV ~ 1 (iV)[iV - 1 + 2(/V - 1) < N > +2(N + 1) < N > 2 ] + • • •} 

G (n e)(0) = 2»N\<N>™ { ATM + 2 < j\f >] 
V 7 (l+2<AT>) Ar + 1 J 

W [TV - 1 + 2iV < TV > +2(/V + 1) < N > 2 } + ■ ■ ■} 

(39) 

p-Gaussian Distribution: 

Per the central limit theorem, the p-Gaussian distribution ( p-normal ) approximation to the 
p-Poisson distribution follows from (18): 

In the p-Poisson distribution function V p (rii,x) for large x = \a\ 2 , we change variables to 
y _ n-fi^ go ag ^ Q measure ^he deviation of n versus the mean p in units of the standard deviation 



p =< a\N\a > 

(40) 

= |«| 2 + |(l-p)-f(l-p) 



12 



a 2 =< a\(N - n) 2 \a > 

(41) 

= \a\ 2 + |(1 - p) 2 + D(l - p)\a\ 2 - \D 2 {1 - p) 2 
Both depend on the positive difference between the n ejD coherent-state-mode probabilities of (14) 

[h_l(\a\ 2 ) - Iz(\a\ 2 )] 
D = DM 2 ) ee P.(|an - P (\a\ 2 ) = + „ > 0- (42) 

For \a\ 2 » 1, 2 



/i ~ |a| 2 + |(1 -p) 



a 2 ~ | a | 2 + i(l-p) 2 



(43) 



since D -> 0. For \a\ 2 -> 0, D -> 1 - + C(|a| 4 ). 

Expanding and using the Stirling approximation, we obtain as in the p — 1 case [8] 

^K^) = ^-^[i-i(f-£) 

-^ + i-*-y 2 {§ + i- 2 ?} + £-S) ( 44 ) 

For V p (n ,y) replace " 72 + f " hi the a -2 coefficient by " — ^ + "• In the bosonic case, 
additional terms through a -4 coefficients are given in [8]. 

As for p = 1, using the first term of this series, we define the p-Gaussian distribution (p-normal) 

G P (n,n,a) ee ^=exp[-i(^) 2 ] (45) 

where /i,a for the pB coherent state are given in (40,41). Unless ^ and/or is large, for \a\ 2 

3 

large the p-Gaussian distribution will be a satisfactory approximation to the p-Poisson when ^- is 
small. 
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